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. Abstract 

I We assume that the Pauli exclusion principle is violated for neutrinos, and cause- 

ys I quently, neutrinos obey the Bose-Einstein statistics. Cosmological and astrophysical 

■ consequences of this assumption are considered. Neutrinos may form cosmological 

Bose condensate which accounts for all (or a part of) the dark matter in the universe. 
i3 I "Wrong" statistics of neutrinos could modify big bang nucleosynthesis, leading to 

>n ■ the effective number of neutrino species smaller than three. Dynamics of the su- 

pernova collapse would be influenced and spectra of the supernova neutrinos may 
change. The presence of neutrino condensate would enhance contributions of the 
Z-bursts to the flux of the UHE cosmic rays and lead to substantial refraction effects 
Ph. for neutrinos from remote sources. The Pauli principle violation for neutrinos can 

<-j ! be tested in the two-neutrino double beta decay. 

> : 

• i-H , 

X 1 Introduction 

What is the next surprise neutrinos will bring us? Pauli has introduced neutrino to resolve 
paradoxes of the beta decay, in particular, an apparent violation of the spin-statistics 
relation. Could the neutrino itself violate this relation? Does the particle invented by 
Pauli respect the Pauli principle? Do we have any indication to that? 

The puzzle of cosmological dark matter (DM) remains with us already for more than 
a half of century but we still do not know what are the constituents of this mysterious 
substance. The commonly accepted point of view is that dark matter is made of new 
elementary particles governed by the laws of the old established physics. Here we will 
explore a different possibility: old particles and new physics. Namely we assume that 
Fermi statistics for neutrinos is violated and, if so, neutrinos with a fraction of eV masses, 
as observed in the oscillation experiments, could condense and make all cosmological dark 
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matter, both cold and (a little) hot. 

We suggest that Pauli exclusion principle is violated for neutrinos and therefore neu- 
trinos obey (at least partially) the Bose-Einstein statistics. Of course, a possibility to 
explain the dark matter is not the only consequence of the spin-statistics violation. One 
can expect some effects of the violation in any environment where large densities or fluxes 
of neutrinos exist. That includes the Early Universe in the epoch of the Big Bang nucle- 
osynthesis, the cores of collapsing stars, etc.. Some consequences could be also seen in 
laboratory experiments (e.g., in the double beta decay). 

Possible violation of the exclusion principle was discussed in a series of theoretical 
papers [I] though no satisfactory model has been proposed so far. (For a critical review 
see ref. j2].) Experimental searches of the Pauli principle violation for electrons [3] and 
nucleons |3] have given negative results. It may happen however that neutrinos due to 
their unique properties are much more sensitive to the violation and it is in the neutrino 
sector the effects can be seen first. Neutrinos may also possess kind of mixed or more 
general statistics than Bose or Fermi ones [S]. 

The assumption of violation of the Pauli exclusion principle reveals immediately a 
number of problems. The spin-statistics theorem follows from the canonical quantization 
to ensure a positive definiteness of energy. It is not clear how to overcome this problem 
and how serious it is for neutrinos. The CPT theorem follows, in particular, from the 
normal relation between spin and statistics, therefore the suggested scenario may also 
violate the CPT theorem. Actually, a possible violation of the latter in neutrino physics 
is under an active study now, see, e.g., refs. 0. Last but not least, the spin-statistics 
violation in the neutrino sector is communicated due to the weak interactions to charged 
leptons and other fermions where the bounds are extremely strong. It is not clear if effects 
considered in this paper are consistent with these bounds, which depends on particular 
mechanism of the violation. 

In what follows we put aside discussion of these problems. Instead, taking pure phe- 
nomenological approach, we concentrate on cosmological and astrophysical consequences 
of the neutrino "bosonization" in an attempt to find interesting observable effects or to 
restrict such a possibility. 

2 Bosonic neutrinos: Context 

The standard electro-weak theory puts the left-handed neutrinos and electrons into the 
same doublet and thus one would expect that neutrinos and electrons obey the same 
statistics. On the other hand, as we know, being the only neutral leptons, the neutrinos 
can have substantially different properties from those of the charged leptons. In particular, 
neutrinos can be the Majorana particles and induce the lepton number violation. The 
difference between the charged leptons and neutrinos is related to breaking of the electro- 
weak (EW) symmetry. The lepton number violation (in the context of seesaw mechanism) 
originates from very high scales. 

Similarly, the neutrino sector might be a source of violation of the spin-statistics rela- 
tion; this can also be connected to EW symmetry breaking and originate from some high 
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mass scale of Nature. One may consider scenario where violation of the Pauli principle 
occurs in a hidden sector of theory related to the Planck scale physics, or strings physics. 
It could be mediated by some singlets of the Standard model - (heavy) neutral fermions 
which mix with neutrinos when the EW symmetry is broken. Since only neutrinos can 
mix with the singlets, effects of the Pauli principle violation would be manifested first 
in neutrinos and then would communicate to other particles. Also one can consider a 
possibility that the messenger of the Pauli principle violation is the light sterile neutrino. 
It has a small mixing with the active components, and this small mixing quantifies the 
degree of violation in the observable sector. In this way a small or partial violation of 
relation between spin and statistics might occur. 

As in the case of lepton number, a violation of the spin-statistic relation for other 
particles can be suppressed by an additional power of a small parameter relevant for the 
violation in the neutrino sector. In fact, the high accuracy of the validity of Fermi statistics 
for electrons [Sj could put a strong bound on a possible "transfer" of wrong statistics from 
neutrinos to electrons. In this connection one can consider a scenario when significant 
effects of the spin-statistics violation develop on the cosmological times or in particular 
environments of the Early Universe. 

A violation of the Pauli principle for neutrinos should show up in the elementary 
processes where identical neutrinos are involved. A realistic process for this test is the 
two-neutrino double beta decay, A — > A' + 29 + 2e~ (or similar with antineutrinos and 
positrons). The probability of the decay as well as the energy spectrum and angular 
distribution of electrons should be affected. 

The difference between neutrino statistics would be strongly pronounced if the prob- 
abilities (see, e.g. [Zj) are proportional to the bi-linear combinations of the type K m K n , 
K m L n , L m L n , where 

K m = [E m — Ei + E el + E u i}~ 1 — [E m — E { + E e2 + E^]" 1 , 
L m = [E m — Ei + E e2 + E ul ] 1 — [E m — Ei + E el + E u2 ] \ (1) 

where E^ is the energy of the initial nuclei, E m is the energy of the intermediate nuclei 
state m, E ei , and E vi are the energies of electrons and neutrinos respectively. The minus 
signs between the two terms in the above expressions are due to the bosonic character of 
neutrinos; in the case of fermionic neutrinos we would have plus signs jZj. For electrons 
we assume the normal Fermi statistics. 

In the case of + — > + transitions the combinations K m and L m can be approximated 

by 

K E e2 — E el + E v2 — E vl E el — E e2 + E v2 — E ul 

(E m -E l + E /2)i ' Lm ^ (E m -E i + E /2y ' [Z) 

whereas for the fermionic neutrinos 

2 

K m w L m m — (3) 

Here Eq/2 = E e + E v is the average energy of the leptonic pair. Appearance of the 
differences of the electron and neutrino energies in @ leads to a suppression of the total 
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probability. It also modifies the energy distributions of electrons. The probabilities of the 
transitions + — > 2 + are proportional to the combinations (K m — L m )(K n — L n ), where 

(K -T \ ~ 2(Ee2 - Eel) (A) 

In the case of the fermionic neutrinos the combination has an additional factor {E v i — 
E„i)/(E m — Ei + Eq/2) and the suppression is stronger. 

A large number of already detected events |H| and especially future measurements 0| 
allow to make precision tests of the Pauli principle. The data seems to exclude the 100% 
violation of Fermi statistics for electronic neutrinos [TUj . Notice however, that relation 
between the statistics of neutrinos and possible (small) violation of the Pauli principle is 
an open issue. 

In what follows we will consider for simplicity an extreme case when neutrinos have 
purely bosonic statistics. In this case the equilibrium density of bosonic neutrinos in phase 
space would be equal to 

f v = (27r) 3 n 5(p) + [exp (E - /i) /T - If 1 , (5) 

where /i is the chemical potential of neutrinos. The first term describes a possible neutrino 
condensate with no being the number density of neutrinos in the condensate. The density 
(jHJ) is somewhat larger than the fermionic one. The thermal part (second term) alone 
gives the energy density of bosonic neutrinos at T ^> m v larger than the fermionic one 
by the factor 8/7. 

Expression (J3J is the equilibrium solution of the kinetic equations since the collision 
integral with such a distribution function vanishes for non-zero n if /x = m v . For smaller 
values of /i, \x < m„, the collision integral vanishes only if n^ = 0. Therefore the condensate 
would be non-vanishing only if the charge asymmetry of neutrinos is so large that the 
maximum chemical potential \i = m u is not sufficient to provide such an asymmetry. We 
will check this explicitly in the next section for the situation during the BBN epoch. 



3 Big bang nucleosynthesis 

The Big bang nucleosynthesis (BBN) is very sensitive to the neutrino statistics. 

Let us first show that the distribution (jSJ) is the equilibrium solution of the kinetic 
equations. For the processes 

e~ + p <-» n + u e , e + + n <-> p + v e (6) 

relevant for the BBN nucleosynthesis, the collision integral is given by the integral over 
the phase space of the function 

F[f] = fj n (l ~ / e )(l - f P ) ~ fef P (l ~ /„)(1 + /„), (7) 

if the invariance with respect to time reversal holds. Here f p , f n and f e are the densities 
of protons, neutrons and electrons correspondingly. We assume that they are described 



4 



by the equilibrium Fermi-Dirac distributions 1 

f F =[exp(E- fJl )/T + l]- 1 , F = e jPl n. (8) 

Notice that the neutrino distribution enters the second them in (J2J) with plus sign. 
Inserting (JHJ) and (JHJ) into (jZJ) we find that F[f] is proportional to 



F oc 



exp^(^ - 1 



1 r 1 

1 - exp — + E v - E e - E p + /i e + ii p - ii n - fi„) 
1 - exp[^(£ n - E e - E p - fi n + /j, e + 



+ (27r) 3 n (5(p) 



(9) 
.(10) 



The first (thermal) term vanishes because of the conservations of energy, E e + E p = 
E n + E u , and the chemical potential: /z e + /i p = /z n + \i v . (The condition for \i is true only 
in equilibrium.) The second term in (fTU|) becomes zero if \i v = m v . Indeed, using the 
conservations of the energy and chemical potential we can rewrite the sum in the exponent 
as E n — E e — E p — fi n + /i e + Hp = —E v + //„ = — m u + /j, u , where the last equality is valid 
since p = 0. So, for fx v = m v , F[f] vanishes and so does the collision integral, even for 
n 7^ 0. Thus, the usual equilibrium distributions remain true despite breakdown of the 
spin-statistics relation. If neutrinos obey mixed statistics their equilibrium distribution 
can be written as a sum of the Bose and Fermi ones at least in the case when a single 
neutrino participates in initial and/or final states of reactions. 

If T-invariance is broken, then a more complicated expression appears in the colli- 
sion integral which still vanishes for the equilibrium distribution because of the S-matrix 
unitarity [TTj. 

Since chemical potential of bosons cannot exceed their mass and neutrino mass is at 
most about eV, the chemical potential of neutrinos would not be essential in the BBN 
range of temperatures, T = (1 — 0.07) MeV. A possible large value of lepton charge 
asymmetry would be hidden in the condensate and would not have a strong impact on 
chemical abundances of light elements. If this is the case, the restrictive BBN bounds O 
18] on the magnitude of the cosmological lepton asymmetry would be inapplicable. 

The equilibrium energy density of the bosonic neutrinos at T ^> m v is 8/7 of the energy 
density of fermionic neutrinos and thus the change of statistics would lead to an increase 
of the effective number of neutrino species at BBN by AN U = 3/7 (for three neutrinos). 
On the other hand, a larger magnitude of the neutrino distribution function and the fact 
that it enters the kinetic equation (see (J7J)) as (1 + /„) instead of (1 — /„) makes the weak 
reactions of neutron-proton transformations (JHJ) faster and the n/p freezing temperature 
becomes lower. The latter dominates and as a result the effective number of massless 
species becomes smaller than 3. 

To estimate the effect we proceed in the following simplified way. The kinetic equation 
which governs the neutron-proton transformation has the form (see, e.g., ref. 



1 Since the violation of the Pauli principle is communicated to other sectors of theory, one may expect 
some deviation from the Fermi-Dirac distribution for other fermions as well, but the effects are probably 
small. 
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where r is the ratio of neutron number density to the total baryon number density, 
r = n/(n + p), x is the cosmological scale factor, = —1.267 is the axial coupling con- 
stant, the function A(T) = B(T) exp (— Am/T), Am = 1.3 MeV is the neutron-proton 
mass difference. The function B(T) is the collision integral containing the equilibrium 
distribution functions of electrons and neutrinos. 

If neutrinos have Bose-Einstein statistics, then the numerical value of B(T) would 
differ from the standard fermionic one. The ratio of fermionic to bosonic -B's at the 
moment of neutron-proton freezing (T « Am/2) is B^/B^ = 0.933. This rise in the 
reaction rate would result in a smaller (n/p)-freezing temperature. The same shift of 
the freezing temperature could be mimicked by a change in the number of the effective 
neutrino species which enters the expression for the Hubble parameter: 

H ~ ^fgl = ^10.75 + 1.75AiV (12) 

One can check that the increase in the reaction rate due to the larger B is equivalent to 
AN V = —0.8. Besides that one should take into account that a real change of the effective 
number of particle species changes not only the neutron-proton freezing temperature but 
also the moment when formation of the light elements begins. It takes place at about 0.07 
MeV and with a smaller number of the particle species this temperature would be reached 
in a longer time and less neutrons would survive the decay. This effect brings about 20% 
into the total impact on AN U during BBN. Since the variation of B does not change the 
time of nucleus formation the total effect of this change on AN U should be —0.6 instead of 
—0.8 found above. Together with the positive contribution from three bosonic neutrinos 
into the total cosmological energy density, which is equivalent to AN U = 3/7, we find 
that the bosonic neutrinos would make the effective number of neutrino species smaller 
by AN V « —0.2. The detailed calculations using properly modified BBN code give [T^j 
AN V = —0.94 for the reaction rate change effect, and 

AN U « -0.51 (13) 

for the total effect. 

The neutrino condensate would not noticeably change the equilibrium value of n/p- 
ratio: 

{n/p) eq = exp [-(Am + /i)/T] (14) 

because the chemical potential of neutrinos is negligible , /i = m v <C T. 

For a positive neutrino asymmetry, i.e., for n„ > n p the condensate can contribute 
to the transition v + n — > p + e, while for a negative asymmetry the contribution of the 
condensate to the (n <-> ^-transformation is kinematically forbidden. These effects are 
automatically included in the case of equilibrium, but when the n/p-ratio is out of equi- 
librium its evolution proceeds somewhat differently from the usual case, and additional 
processes with condensate should be taken into account. This would make AN U a few per 
cent more negative. 

The negative value of AN U found above seems to be in concordance with the data on 
the light element abundances. According to the recent analysis ^H] 2.67 < N u < 3.85 
at 68% CL, while ref. [Hj presents the result AA^ = — 0.37lo!i?- Anyhow, taken at the 
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face value, the observed abundances of 4 He and 2 H seem to be in contradiction with 
the standard BBN calculations, especially if lower values of the observed abundances are 
taken. At the same time AN V = —0.2 opens some room for additional light degrees of 
freedom. 

If neutrinos have mixed statistics, e.g., the Pauli principle is broken for sterile neutrinos 
which mix weakly with the standard ones, the effect of "wrong" statistics on N v would 
be further diminished. 



4 Bosonic neutrinos and dark matter 

If neutrinos obey the Bose-Einstein statistics, the Gunn-Tremaine JH] lower bound on 
their mass would be inapplicable. This bound is based on the Fermi statistics which does 
not permit to have too many neutrinos in a galaxy. Thus, to make all galactic dark matter, 
neutrinos (or any other fermions) must be sufficiently heavy, m„ ~ 100 eV. However, this 
large mass contradicts the Gerstein-Zeldovich upper bound JH]- Bosonic neutrinos can 
be arbitrary light, even as light as axions, but still make all cosmological dark matter. 

To play a noticeable role in the cosmological large scale structure (LSS) formation 
and, possibly, to make all the observed cold dark matter the neutrino condensate should 
be sufficiently large. For m u rs 0.1 eV the required energy density Q c = 0.25 or p c = 
1.25 keVcm -3 can be achieved if the number density of neutrinos in the condensate is 
1.25 • 10 4 cm -3 . It is approximately 100 times larger than the number density of neutrinos 
plus antineutrinos for any single neutrino species from the cosmological thermal bath. To 
create such a large number of cold neutrinos the cosmological lepton asymmetry must be 
much larger than 1. 

The large lepton asymmetry can be created in a version of the Affleck-Dine sce- 
nario [201 • Suppose there is a scalar field, x> with non-zero lepton charge. This field 
can acquire a large vacuum expectation value during inflation. The potential of x self- 
interaction, U(x), is supposed to break the leptonic charge conservation, and therefore 
U(x) is not invariant with respect to phase transformation, x ~^ X ex P(^). Usually it is 
assumed to have the form 

U( X ) = ml\ X \ 2 + mix 2 + mf X * 2 + A^ 4 + W + Kx*\ (15) 

where mi and Ai are real but m 2 and A 2 may be complex so that the C-invariance is 
broken. 

If the mass parameters m 12 are small in comparison with the Hubble constant during 
inflation, Hj, then x 2 could acquire a large vacuum expectation value 

if 4 H 2 

(X 2 )~^ ox~-L. (16) 
mf Ai 

Moreover, along possible flat directions of the potential (valleys) x 2 would rise with time 
as Hft (see the discussion and the list of references, e.g., in [21] ) . 
The initial leptonic charge density of x is equal to 

L x = e\x\ 2 , (17) 
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where 9 is the phase of the field, x — \x\ exp(z0). The energy density of x is about the 
kinetic energy, # 2 |x| 2 > and the potential energy is normally of the same order of magnitude. 

When inflation is over and the Hubble parameter becomes smaller than the mass of 
X, the field starts to relax down to zero and in the process of relaxation it transmits the 
accumulated leptonic charge into leptonic charge of the decay products. We assume that 
X decays into two neutrinos, x ~* vv i (arid not antineutrinos) , i.e., the leptonic charge of 
X is two. In this case the lepton number density of the produced neutrinos would be 

n L = L x . (18) 

If the decay and the subsequent thermalization are fast, the energy density of the created 
plasma becomes 

PT = e 2 \ X \ 2 = (n 2 gj30)T\ (19) 

where g* is the number of relativistic degrees of freedom in the plasma. 

It follows from (|17j) . ()18j) and that the charge asymmetry of the plasma could be 

/3 L = n L /T 3 ~ ^JxTo- (20) 
For the initial values 9 ~ Hi and x ~ H 2 /m, the asymmetry equals 

Pl ~ -, (21) 
m 

and it does not look unreasonable, e.g., that (3l > 100. 

After x decayed into neutrino pairs, the plasma thermalized through the reactions 
vv — > vvll, etc.. Complete thermalization would produce the equilibrium distribution 
of neutrinos given by eq. (J5)l. Because of a low rate of reactions with zero momentum 
neutrinos the complete equilibrium with condensation at p u = is never reached, still, 
according to the estimates presented below, neutrinos accumulate sufficiently close to 
p = and the condensate practically forms. 

The kinetic equilibrium could be achieved through elastic scattering reactions v + 
I — > v + /. If the initial temperature was above the electroweak scale and therefore the 
intermediate W and Z bosons could be considered massless, the low energy band of the 
spectrum of neutrinos would be filled in with the rate 

T ~ g 4 E u ~ a 2 E u , (22) 

where g is the electroweak coupling constant, E u is the neutrino energy and a = 1/137. 
At temperatures below the masses of W and Z bosons, the rate becomes much lower: 

T ~ G 2 F T 4 E U , (23) 

where Gp is the Fermi coupling constant. 

Requiring r > H, we conclude that neutrinos "condense" down to energies E u /T ~ 
lOT/mpi if the process occurs before the EW phase transition, i.e. if x decays at tem- 
peratures T > mwz- When T drops below mw,z the cooling would become less efficient. 
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Effective energies to which neutrinos could cool down in the condensation process in- 
creases up to I {T 2 m P i). Though complete condensation with p = could probably 
never be achieved, an accumulation of neutrinos at low energy part of the spectrum seems 
to be efficient enough to create cosmo logical cold dark matter of neutrinos 2 . 

In contrast to the case of large charge asymmetry carried by fermions (when the 
spectrum is populated by energetic particles), in the bosonic case the charge asymmetry 
is stored in the condensate, while chemical potential \x = m u is very small and the energy 
density is dominated by relativistic particles with practically vanishing charge asymmetry. 
It remains true till temperature dropped down below m v . After that the energy density 
would be dominated by non-relativistic condensate with large leptonic charge. 

The energy density in the neutrino condensate equals 

Pcond = m„n . (24) 

Requiring that at the BBN epoch p con d < O.lp^, where pi is the energy density in the 
thermal bath, we find the bound 

n < 0.3 n ~ 10V. (25) 

m v 

Here is the neutrino number density in the thermal bath. 

If the equilibrium concentration (J3J) is not achieved and the condensate is not formed 
completely, so that p ^ 0, the BBN gives the following bound: 

T 

fl 

|p| < 0.3 Tbbn ■ (26) 
n 

For nl/riQ ~ 10~ 2 required to explain the dark matter we find |p| < 3 ■ 10 3 eV. 



The large scale structure formation with cold dark matter composed of the neutrino 
condensate would well fit the observed picture. The analysis of structure formation with 
the Bose condensate of normal bosons was done in ref. [22] and it can be extended to the 
case of neutrino without significant modifications. 

One can be less ambitious and admit two coexisting forms of the cold dark matter: 
e.g., the usually accepted lightest supersymmetric particle (LSP) or axions, and the Bose- 
condensed neutrinos. It would make the scenario considered here less vulnerable but 
simultaneously less predictive. 

Notice that because of smallness of the neutrino chemical potential a large lepton 
asymmetry in primordial plasma cannot be transformed into the baryon one by elec- 
troweak processes. Indeed, the chemical potential of neutrinos is tiny and in equilibrium 
the same must be true for chemical potentials of quarks. Hence the baryon asymmetry 
generated by sphalerons in equilibrium could be at most of the order m v /T ~ 10~ 12 . 

2 In the case that sterile neutrinos condense, the necessary for DM magnitude of the condensate should 
be scaled by the ratio of sterile/active masses. 
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5 Astrophysical consequences 



The neutrino statistics plays the key role in the environments where neutrinos form dense 
degenerate gases. 

Direct test of the "bosonic" nature of neutrinos can be probed by precise measurements 
of the neutrino energy spectrum from supernova. Instead of the Fermi-Dirac spectrum 
with pinching the distribution would be the Bose-Einstein one with some pinching effect 
too. So, generically, the spectrum of bosonic neutrinos should be narrower. To establish 
the difference one needs to measure the spectrum both in the low, E < 3T, and in the 
high, E > 3T energy parts. Also pinching effect should be quantified rather precisely. 

Violation of the Pauli principle can influence dynamics of the SN collapse. According 
to the usual scenario in the initial stages (formation of the hot proto-neutron star) the 
neutronization leads to production of high concentration of the electron neutrinos which 
are trapped in the core. The chemical potential of these neutrinos (due to the Pauli 
principle) can reach 70 - 100 MeV. These neutrinos heat the medium and diffuse from 
the core. Violation of the Pauli principle allows for the neutronization neutrinos to be 
produced with lower energies. These neutrinos escape easier the star leading to faster 
cooling and lower central temperatures. Also the evolution of the lepton number would 
change. 

The presence of the neutrino condensate with large lepton number in the Universe 
may have a number of observable consequences. 

High neutrino density in the condensate (especially if an additional clustering occurs) 
enhances rate of the Z°-bursts produced by annihilation of the ultra high energy (UHE) 
cosmic neutrinos on the relic neutrinos [2Hll2l]- This in turn, enhances production of the 
UHE cosmic rays, and may help to explain the cosmic ray evens above the GZK cut-off. 

The asymmetric neutrino condensate may produce strong refraction of the high energy 
neutrinos from remote sources (active galactic nuclei, gamma ray bursters). Apart from 
lensing, one may expect a substantial impact on neutrino oscillations [T2] . 

Since the density of dark matter in galaxies is about 6 orders of magnitude larger than 
their average cosmological energy density, a condensation of cold neutrinos around the 
Earth might have an effect on the end point of the beta decay spectra, in particular, in 
the tritium decay experiments on search for neutrino mass. 



6 Conclusions 

If the Pauli principle is violated in Nature, neutrinos may be the first messengers of that. 
Indeed, neutrino properties may be related to physics at high energy scales and only neu- 
trinos (due to their neutrality) can mix with particles of the hidden sector (singlets of the 
standard model symmetry group). 

We have considered the cosmological and astrophysical consequences of the Pauli princi- 
ple violation assuming that neutrinos obey the Bose-Einstein statistics. Such neutrinos 
would have the Bose-Einstein energy distributions and might form a condensate. 
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We have found that in the case of bosonic neutrinos, the effective number of the neutrino 
species during the BBN epoch is reduced by about AN U = —0.5, thus slightly facilitating 
an existence of new degrees of freedom, e.g., sterile neutrinos. 

If in the Early Universe, a large lepton asymmetry was created, bosonic neutrinos should 
condense in the process of the cosmological cooling. Because of that the BBN bound on 
the neutrino density can be avoided, and so a large neutrino concentration and energy 
density at a later epoch become possible. In particular, it allows to explain all or a part 
of the cosmological dark matter. It opens a possibility of large lepton asymmetry of the 
Universe. 

Bosonic neutrinos would influence dynamics of the stellar collapse and lead to modifica- 
tion of the energy spectra of SN neutrinos. 

A large concentration of relic background neutrinos in the condensate at the present epoch 
would enhance the rate of the Z°— bursts. It can produce strong refraction effect on neu- 
trinos from remote sources, modifying the oscillation pattern. 

A study of the double beta decay, in particular, of the energy spectrum and angular 
distribution of electrons from the decay can provide sensitive test of the Pauli principle 
violation and statistics of neutrinos. 

Note added: After submission of the first version of our paper, we became aware of the 
paper [23] in which a possibility of bosonic statistics for neutrinos has been considered and 
its effects on the BBN have been studied. According to the change of neutrino statis- 
tics is equivalent to the decrease of number of the effective neutrino species AN V = —0.74. 
This conclusion agrees qualitatively with our results of sec. 3, though quantitatively we 
find smaller decrease. 
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